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We study the QCD vacuum structure under the influence of an elecfromagnefic field wifh a nonzero 
second Lorentz invarianf I 2 = E ■ B. We show fhaf fhe presence of I 2 can induce neufral pion (tt^) 
condensafion in fhe QCD vacuum fhrough the electromagnetic triangle anomaly. Within the frame¬ 
works of chiral perturbafion fheory at leading small-momenta expansion as well as the Nambu-Jona- 
Lasinio model at leading 1/Nc expansion, we quantify fhe dependence of fhe condensate on 12- 
The stability of fhe 7r*^-condensed vacuum againsf fhe Schwinger charged pair producfion due fo 
elecfric field is also discussed. 

FACS numbers: 12.38.Aw, 12.39.Fe, 11.30.Rd 


INTRODUCTION 


The vacuum structure of Quantum Chromodynamics 
(QCD) has been intensively studied over many years 
and a variety of interesting properties have been ex¬ 
tracted. For example, it is well-known that the normal 
QCD vacuum possesses the spontaneous (approximate) 
chiral symmetry breaking (CSB) which can be character¬ 
ized by the scalar quark-antiquark condensate {^tp) 7^ 
0; while the pseudoscalar condensates like {ipij^Tstp) 
with T3 the third Pauli matrix are not permitted — an 
assertion known as the Vafa-Witten theorem [1]. 

Under proper conditions, the QCD vacuum structure 
can be altered, usually through phase transitions. Exam¬ 
ples include the pion superfluid phase at high isospin 
chemical potential [2-4] and the color superconducting 
phase at high baryon chemical potential [5-7]; in the 
first case the isospin symmetry is broken while the chi¬ 
ral symmetry is restored and in the second case the color 
symmetry is broken while the chiral symmetry may or 
may not be broken. 

Strong electromagnetic (EM) fields provide another 
way to modify the QCD vacuum. Strong magnetic 
fields may exist in compact stars [8-10], heavy-ion colli¬ 
sions [11-13], and early universe [14-16]. Strong electric 
fields may also be accessible in heavy-ion collisions [17- 
19]. The presence of the magnetic field is known to en¬ 
hance the scalar {^tp) condensate at zero quark chemical 
potentials which is called the magnetic catalysis of the 
CSB [20-23], while the effect of the magnetic field on the 
CSB at finite temperature or chemical potentials shows 
novel features which are still not fully understood [24- 
29]. It was also proposed that when \eB\ > nip with 
nip the mass of the rho meson, the QCD vacuum will 
become a superconductor due to charged rho condensa¬ 
tion [30]. The effect of electric field on the CSB was also 
studied [31-33] and it was found that the electric field 
always tends to break the scalar quark-antiquark pairs 
and thus weaken the CSB. 

In this Letter, we study the combined effect of the elec- 
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FIG. 1: The electromagnetic triangle anomaly. 


trie and magnetic fields on QCD vacuum. We focus on 
the EM-field configuration where the second Lorentz in¬ 
variant I 2 = E ■ B is nonzero. We show that the pres¬ 
ence of I 2 provides a parity-odd environment in which 
the otherwise-forbidden (by Vafa-Witten theorem) neu¬ 
tral pion condensation can occur via the EM triangle 
anomaly (see Eig. 1). We calculate the condensate 
within both the chiral perturbation theory (ChPT) and 
the Nambu-Jona-Lasinio (NJL) model frameworks and 
identical result is obtained. 
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Nc 


^ ' Sgn(t2) 


for|J2|<t^, 
for |t2| > t2' 


( 1 ) 


where m* = {n^)^ + cr^ with a ~ {^tp) and ~ 

{{pij^T^ip) the scalar and neutral pion condensates in the 
presence of the EM field, qu = 2e/3 and = —e/3 
are charges of u and d quarks with e > 0 the proton 
charge, and I 2 = 4:71^f^nt^/[Nc{qu ~ ‘?d)] is ^ critical 
value for |t2| above which the condensate ni* is wholly 
contributed by mode. Equation (1) shows that the 
I 2 rotates the chiral condensate from the u-direction 
to the 7r*^-direction with the rotation angle given by 

(p = sin~^{ tz^/ m*) if |t2| < ‘P ~ if 

\h\ > I 2 (see Eig. 2). 

Throughout this article, we restrict to the zero quark 
chemical potential and zero temperature case. Before 
we proceed, we note that the neutral pion condensa¬ 
tion we will study is different from the disoriented 
chiral condensate (DCC) [34] proposed to possibly oc¬ 
cur in heavy-ion collisions, as the latter is a far-from- 
equilibrium, transient, phenomenon while the former is 
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FIG. 2: Illustration of the chiral rotation due to the neutral pion 
condensation induced by EM triangle anomaly. 


CHIRAL PERTURBATION THEORY CALCULATION 

We start with the two-flavor ChPT described by the 
Lagrangian 


£ — £o + £wzw/ 


(3) 


where Cq is the usual chiral Lagrangian given by (we 
keep only O(p^) terms) 


/■2 

£o = ^tr 


D^U+Di‘Lf- 


m\{U + U^) 


(4) 


and the Wess-Zumino-Witten term C^zw is given 
by [41-45] 


a nearly equilibrium, (quasi-)stationary, state. We point 
out that, in the context of DCC, it was also found that the 
large EM fields in heavy-ion collisions could also possi¬ 
bly catalyze the formation of the DCC [35, 36]. 


INVALIDATION OE THE VAEA-WITTEN THEOREM 

The EM-triangle-anomaly induced neutral pion con¬ 
densation does not contradict the Vafa-Witten theorem. 
The validity of the Vafa-Witten theorem relies on the 
positivity, detH > 0, where D is the Euclidean Dirac 
operator [1,37]. However, this is not the case in the pres¬ 
ence of a real electric field. To see this, we write down T) 
explicitely. 




£wzw = (QL^,L«L^ + QRyRaRfi) 

-iFocfiTv], (5) 

q2(L, + Rv) + i [qUQU^L, + QU^QUR 


Tv = tr 


If, = Ud^U\ 
Rf, = df,U^U. 


( 6 ) 


In the above, the covariant derivative is given by 

Df,U = d^U + Af,[Q,U], (7) 


and Lf is a 2 X 2 unitary matrix representing the chiral 
fields, for which we choose the Weinberg parametriza- 
tion. 


T’ — + Q74A0 + iQ'YiAi + M, ( 2 ) 

where (// = 1 — 4) is the gluon field, (a = 0 — 3) 
is the EM field, and M is the mass matrix. We choose 
the gauge such that E = — VAq and A, is time in¬ 
dependent. Eor two-flavor case, the charge matrix is 
Q/e = diag(2/3, —1/3) = 1/6-|-T 3/2 and the mass 
matrix is M = mq. The crucial observation is that Q is 
not traceless and thus its role is similar with the baryon 
chemical potential }i. It is well-known that the presence 
of J.I destroys the positivity of detH, and so does Aq. In 
fact, when Aq 7 ^ 0 there does not exist a matrix T such 
that T~^'DT = to guarantee the positivity. Choosing 
another gauge, e.g. Aq = 0 , does not change the conclu¬ 
sion. Because in this case. A, must depend on time, e.g. 
for constant E, A, = — E,f, which after Wick rotation 
spoils again the positivity [38]. Note that if qu = —qd, 
deiV is semi-positive [39]. 

The invalidation of the Vafa-Witten theorem in the 
electric field has important physical contents. Physi¬ 
cally, in the presence of E, the Dirac determinant has to 
be complex or at most semi-positive because it always 
receives contributions from the Schwinger poles which 
represent the particle-antiparticle pair creation due to 
E [40]. 


U=^{s + iT-t), ( 8 ) 

Jn 

where the fields s and t fulfill the constraint = fn- 

Now we consider constant EM field and make the ap¬ 
proximation such that we treat all the chiral fields in 
(3) as uniform condensates. We note that the a and rfi 
in Eq. (1) are proportional to s and respectively. In 
the absence of the EM field, it is well-known that the 
CSB takes place along the n-direction so that s = fn 
and t = 0. The EM field with a configuration where 
I 2 = E ■ B ^ 0 is coupled to the rfi field via the trian¬ 
gle anomaly and raises the possibility of condensa¬ 
tion. If this happens, the CSB will be rotated from the 
n-direction to Ti'^-direction. To explore this possibility, 
we make the following ansatz for the condensates, 

s = fn COS (p, ts = fn sin (p, (9) 

h = t2 = 0, (10) 

where without loss of generality, we choose cp € 
[—7t/ 2, 7r/2] so that s > 0. The Lagrangian £ then be¬ 
comes 

C{(p) = f^ml cos (p + {Q%)(p, (11) 
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where we have omitted a total derivative term. As ob¬ 
viously, the second term is responsible for the triangle 
anomaly effects such as rfi —> 27 . It is easy to minimize 
— C{(p) and obtain 


sin^ 


‘P 


Nch 


Nch 




sgn {h) for 


2 ' 


h(Q^T3) 

iql - 

\h\ 


for |f 2 | < I 2 , 

>I2, 


( 12 ) 

(13) 


where = 47 T^/^m^/[Nc(g^ — q^)]. This is equivalent 
to Eq. (1). 

Several comments are in order. (1) In the chiral limit 
mjr —> 0 , thus I 2 —t 0 , an infinitesimal I 2 will drive 
the maximum condensation, namely, the CSB is com¬ 
pletely driven to tt*^- direction once I 2 is applied. (2) The 
neutral pion condensation induced by I 2 is merely a chi¬ 
ral rotation; it is not a phase transition as no symmetry 
is broken along this process. (3) As a low-energy effec¬ 
tive theory of QCD, the ChPT is reliable only when all 
the physical parameters, rtiji, p, ^y\eE\, are much 

smaller than the typical hadronic scale ^ 1 GeV. We 
can estimate that 






—— 1 


(14) 


which justifies the validity of our result for I 2 not too 
larger than (4) As well-known, when \eE\ ^ 

(but not larger than A^), the vacuum will be unstable 
because of the Schwinger pair production of 71+ n~. We 
will discuss this issue in another section. 


where ip = {u,dY represents the two-flavor quark 
fields, mg is the current mass of quarks, G is the four- 
fermion coupling constant, and 

-|- zQA^, (16) 

is the covariant derivative. Note that, in addition to the 
coupling constant G, the NJL model contains implicitly 
another parameter, the ultraviolet cutoff A, which spec¬ 
ifies the applicable region of the model. The values of 
G and A are determined by fitting NJL predictions with 
the physical hadronic observables like fn and in 
the normal vacuum [46]. 

The NJL model was widely used to study a number 
of non-perturbative properties of QCD, especially those 
related to the chiral symmetry. In fact, the NJL model 
shares the same global symmetries with QCD. At the 
mean-field level where only the leading-order terms in 
the 1 / Nc expansion are kept. 


>Cnjl = -m- i7r°75T3)i/^ - ^ , (17) 

where m = mQ + cr is the constitute mass of quarks and 
a = —2G{xpq}) and = — 2 G( 1 /^ 775 X 3 i/;) represent the 
sigma and neutral pion condensates. The quark fields 
can then be integrated out analytically, which is most 
conveniently done in Euclidean path integral and the re¬ 
sultant thermodynamic potential reads 


O = 


(m 


nio)^ + ( 7 t °)2 

4G 


:^TrlnS-\ 

1/4 


(18) 


NAMBU-JONA-LASINIO MODEL CALCULATION 

The ChPT is built on hadronic degrees of freedom, we 
now use the NJL model which is based on quark degrees 
of freedom to test the consequences of the ChPT. The 
Lagrangian of the two-flavor NJL model is given by [46] 

>Cnjl = + G[(#)^ -h (Tpi-YsTipf], (15) 


where V 4 is fhe Euclidean spacetime volume, Tr acts on 
all physical spaces, and the quark propagator is 

S(x,x') = —(z0 — m — 771 *^ 75 X 3 )“^^^^^ (x — x'). (19) 

It is diagonal in flavor space which allows us to calculate 
it analytically by adopting the Schwinger proper-time 
formalism [40]. The result is 


_ 1 pCQ px A j 1 

Sf{x,x') = [-+"^-sgn(<?/);V 7 T°] 


X exp I — {[np- + (7r‘^)^]s -|- ^{x — x')qfFcoth{qj:Fs){x — x') -|- ^q^crFs'^ 


-{q/sfh 


Imcosh [iqfs{Ii -|- 2il2p^^)' 


where f = u,d, Ii = — E^ is the first Lorentz in¬ 


variant, and we have suppressed some of the Lorentz 
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scripts, e.g., aF should be understood as cr^^'^Fy^ with 
The true vacuum corresponds to the 
global minimum of O. In the case that = 0 and 
I 2 > 0, we have numerically minimized H and obtained 
the corresponding m and n^, see Fig. 3. Once I 2 is turned 
on, finite condensate arises and sigma condensate is 
suppressed — a feature representing a chiral rotation 
from n-direction to 7r*^-direction; when ^ 0.28 GeV, 


the CSB is driven by condensate solely. In addition, 
the total condensate m* = is weakened 

and finally killed by I 2 due to the presence of the elec¬ 
tric field [32,33]. 


To gain more insights, we write down the gap equa¬ 
tions which is valid when 1 121 < ^2 below): 


m — niQ 
IG 

2G 


mNc /■“ _[,„2+(„0)2]g qjhRecoshlqfSy'h+lih] j^o ^2 2^r 

^71^ fk,dJo lmcosh[qfSy/h+2il2\ 

n°Nc ^ r, _[„ 2 +(^ 0 ) 2 ], ^?J/ 2 Recosh[^/ysV-fi + 2 /I 2 ] m ,2 2^r 

Imcosh[qysVfi+2fl2] 4712 m 2 + (710)2^'^“ 


( 20 ) 

( 21 ) 


By eliminating m, we analytically obtain 


STABILITY OF THE 7r“-CONDENSED VACUUM 


770 


Nc 2G, 2 


7Id)h- 


( 22 ) 


By assuming the NJL version of the Gell-Marm-Oakes- 
Rermer relation, tn^f^ = mom* {2G)~^ [46, 47], Eq. (22) 
immediately reduces to Eq. (1) for IJ 2 I < which is 
independent of the model parameters G and A. 

Two comments: (1) If mg = 0 or m„ = 0, the above 
formula is divergent. This corresponds to the chiral limit 
case where once I 2 is turned on the CSB is immediately 
rotated to the 7r*^-direction. (2) When |f 2 | > 
the gap equations is not applicable because m is already 
zero. This situation is consistent with the ChPT result. 
In both NJL and ChPT frameworks, when 1121 > I 2 , the 
CSB is found to be totally driven by condensation. 
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FIG. 3: The constitute quark mass m and condensate as 
functions of obtained by minimizing the thermodynamic 
potential Q at Ij = 0 in NJL model. The presence of I 2 always 
tends to diminish n condensate while drive condensation. 


The neutral pion condensation found above requires 
a nonzero electric field. As we know, strong electric field 
can efficiently generate particle-antiparticle pairs (in the 
confined phase, mostly 71+ and n~; in the deconfined 
phase, quarks and antiquarks) through the Schwinger 
mechanism and drive the vacuum unstable. So it is nec¬ 
essary to analyze that rmder which condition the tt*^- 
condensed vacuum is stable. For this purpose, we need 
to calculate the pair creation rate. This can be achieved 
within the ChPT framework by using the Schwinger 
proper-time method, see, e.g.. Refs. [33, 48]. We will not 
repeat their calculations here, but just quote their results 
relevant to our discussion. 

In case that I 2 yf 0, we can always find a frame in 
which E and B are parallel. We will work in such a 
frame. The imaginary part of the leading order effective 
Lagrangian, which is responsible for the 7r+ n~ produc¬ 
tion, reads 


Im£ = 


e^lEBl 


167r2 


~ (-l)«+i 

L—7— 


n,l—l 


—nn 


m2+|eB|(2/-l) 

leEl 


(23) 


where n runs over all Schwinger poles and I runs over 
all Landau levels. The first term in the summation over 
n defines the tt+tt” production rate per unit volume 
which reads, 

e^lEBl ~ 

It is clear that the presence of the magnetic field effec¬ 
tively enhances the mass of charged pions. 


— n 


mi, + \eB\{2l - 1) 


eE 


.(24) 


= ml+\eB\, 


(25) 
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and therefore as long as 

|eE| <C + |eB|, (26) 

the rate of tC production is strongly suppressed. In 
this case, the pair-producing instability takes place over 
a very long time f ~ 1 / {w„+„- V) {V the volume) and 
it is meaningful to consider the "static" property of the 
vacuum. The "equilibrium state" neutral pion conden¬ 
sate we have studied so far is sensible when the condi¬ 
tion (26) is satisfied. 


SUMMARY 

In this article, we propose a mean to tune the 
QCD vacuum by applying parallel electric and mag¬ 
netic fields. The underlying mechanism is the electro¬ 
magnetic triangle anomaly, through which the second 
Lorentz invariant I 2 = E ■ B of the EM field can induce 
a neutral pion condensate in QCD vacuum and lead to a 
chiral rotation from the isosinglet scalar cr-direction to¬ 
ward the isotriplet pseudoscalar 7 T*^-direction. By adopt¬ 
ing the chiral perturbation theory at leading order plus 
the Wess-Zumino-Witten term and the Nambu-Jona- 
Lasinio model at the mean-field approximation (leading 
order in 1/Nc expansion), we show that the sine of the 
rotation angle is universally given by Eq. (l).Eor weak 
electric field, our results may be testable in lattice QCD 
by adopting similar strategies as the ones used in lattice 
simulations at finite baryon chemical potential [49]. 

Our finding may have experimental implications in 
heavy-ion collisions. In ultra-peripheral heavy-ion col¬ 
lisions the hot quark-gluon matter may not form but 
quite large I 2 with opposite signs above and below the 
reaction plane may exist [12] which can induce transient 
71^ condensation. The diphoton spectrum may be used 
to distinguish the 7r*^-condensed vacuum from the nor¬ 
mal vacuum, because in the Tr’^-condensed vacuum the 
a and excitations are mixed which will change the 
diphoton emission rate and spectrum. We will study 
these experimental signatures in near future. 
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